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ABSTRACT 


The process of cold flat plane-strain rolling of 
rigid-plastic strain hardening materials is analyzed by Galerkin 
finite element technique coupled with the direct penalty method at 
steady state. Analysis is carried out both without and with strain 
hardening and the analytical results for normal pressure 
distribution, roll force and roll torque with strain hardening are 
compared with the experimental results. A detailed parametric 
study is carried out to show the effect of the process variables 
viz. reduction ratio, R/h^ and the coefficient of friction on the 
important rolling design parameters of normal pressure 
distribution, roll separating force and roll torque. The 
percentage under estimation of these rolling parameters with the 
assumption of perfectly plastic material is highlighted. A study 
of contours of normal strain -rate, equivalent strain -rate, 
deviatoric stress and equivalent strain for four selected sets of 
process variables is presented. It is concluded that the present 
model is capable of estimating the important rolling parameters 
for a range of process variables and easily accommodates 
experimentally measured variable coefficient of friction for- 
estimating more accurate results. 


iii 



ACKNOULEDGEHENTS 


I acknowledge with aincerity and deep aense of gratitude the 
expert guidance and continuoua encouragement provided by 
Dr. G.K.Lal and Dr. P.M. Dixit throughout the courae of thia 
theaia. I alao thank them for their contr ibut iona made by the way 
of advice and criticiam at varioua atagea of thia work. 

I thank the entire ataff of the computer center, I IT Kanpur, 
for enabling me to carry out the computationa without any major 
interruption, albeit aome occaaional ahutdowna. 

I wiah to thank all my frienda, neighbor and acquaintances in 
the campus for making my stay at I IT Kanpur memorable. 


Regalia Srinivasa Prakash 

April, 1992 
I . I . T . Kanpur . 


iv 



TO 

MY BELOVED PARENTS 



CONTENTS 


Page No. 

List o£ Tables vill 

List of Figures 

List of Notations xiii 

CHAPTER 1 INTRODUCTION 

1.1 General Introduction 1 

1.2 Mechanics of Rolling Process 1 

1.3 Previous Work 3 

1.4 Objectives of the Present Work 6 

1.5 Scope of the Present Uork 8 

1.6 Plan of the Thesis 10 

CHAPTER 2 MATHEMATICAL MODELLING OF PLANE-STRAIN COLD 

ROLLING PROCESS 

2.1 Introduction 11 

2.2 Governing Equations 11 

2.3 Equations for Steady State 

Two Dimensional Problem 13 

2.4 Boundary Equations 15 

2.5 Non-Dlmenslonalization 17 

2.6 Galerkin Finite Element Formulation 19 

2.6.1 Galerkin Ueak Formulation with 

Direct Penalty Method 19 

2.6.2 Finite Element Approximation 22 

2.7 Finite Element Equations in Local Variables 25 

Vi 



2.8 Glob&l Aaaeflibly 28 

2.9 Application o£ Boundary Conditiona 30 

2.10 Formulation for Strain Hardening 33 

2.11 Implementation 36 

2.12 Evaluation of the Secondary 
CHAPTER 3 RESULTS AND DISCUSSION 

3.1 Validation 48 

3.2 Parametric Study 49 

3.2.1 Effect of Reduction Ratio and R/h^ 49 

3.2.2 Effect of Coefficient of Friction 51 

3.3 Effect of Strain Hardening 52 

3.4 Study of Contours 54 

3.5 Conclusions 57 

3.6 Scope for Further Research 57 

REFERENCES 84 

BRIEF BIBLIOGRAPHY 87 


vii 



LIST OF TABLES 


page no 


Table 

3 . 1 

79 

Table 

3.2 

79 


vlli 



LIST OF FIGURES 


Figure No. 
Fig. 2.1 
Fig. 2.2 
Fig. 2.3 

Fig. 2.4 
Fig. 2.5 
Fig. 2.6 
Fig. 2.7 
Fig- 3.1 

Fig. 3.2 

Fig. 3.3 

Fig. 3.4 

Fig. 3.5 

Fig. 3.6 


Title Page No. 

The Flat Rolling Process 

The Eulerian Frame 43 

The Boundary of the Domain and the 44 


Boundary Conditions 
A Typical Finite Element 
A Typical Boundary Element 
The Boundary Tractions 
The Finite Element Nesh 

Comparison of Analytical and Experimental 59 

Roll Separating Force for St eel (o'^ = 324 MPa) 
Comparison of Analytical and Experimental 59 

Roll Separating Force for St eel (o'^ = 324 MPa) 
Comparison of Analytical and Experimental 60 


Roll Separating Force for Aluminium at 
R/h^ = 12.5 

Comparison Analytical and Experimental 60 

Roll Separating Force for Aluminium at 
R/h^ = 39 

Comparison of Analytical and Experimental 
Roll Separating Force for Copper at 
R/h^ = 12.5 

Comparison of Analytical and Experimental 51 

Roll Separating Force for Copper at 


ix 



50 


= 

Fig. 3.7 Comparison of Analytical and Experimental ^ 

Roll Torque for Steel = 358 MPa ) 

Fig. 3.8 Comparison of Analytical and Experimental 52 

Roll Torque for Steel = 324 MPa ) 

Fig. 3.9 Comparison of Analytical and Experimental 63 

Roll Torque for Aluminium at R/h^^ = 12.5 
Fig. 3.10 Comparison of Analytical and Experimental 63 

Roll Torque for Aluminium at R/h^ = 39 
Fig. 3.11 Comparison of Analytical and Experimental 

Roll Torque for Copper at R/h^ = 12.5 
Fig. 3.12 Comparison of Analytical and Experimental 64 

Roll Torque for Copper at R/h^ = 50 
Fig. 3.13 Comparison of Analytical and Experimental 65 

Normal Pressure Distribution for 
Aluminium at R/hj^ = 12.5; %r = 14.17; f = 0.10 
Fig. 3.14 Comparison of Analytical and Experimental 65 

Normal Pressure Distribution for 
Aluminium at R/h^=39; %r=33.75; f=0.10 
Fig. 3.15 Comparison of Analytical and Experimental 66 

Normal Pressure Distribution for 
Copper at R/hj^ = 12.5; %r=l4-.20; f = 0.06 
Fig. 3.16 Comparison of Analytical and Experimental 56 

Normal Pressure Distribution for 
Copper at R/hj^ = 50; %r= 2T .76 ; f = 0.06 

Fig. 3.17 Normal Pressure Distribution Along the 57 

Arc of Contact at Different Reductions 


X 



Fig. 3.18 


68 


for Steel at R/h^= 65 ; f = 0.06 

Normal Preaaure Diatrlbution Along the 

Arc of Contact at Different Valuee of 
R/hj^ for Steel at kc= IQ ; f=0.06 

Fig. 3.19 Shear Streaa Diatrlbution Along the 69 

Arc of Contact at Different Reductions 

for Steel at R/h^= ( K ; kr= \6 ; f = 0.06 

Fig. 3.20 Shear Stress Distribution Along the 70 

Arc of Contact at Different Values of 
R/h^ for Steel at %r= I 6 ; f=0.06 
Fig. 3.21 Variation of Roll Separating Force 71 

with R/hj^ at Different Reductions 
for Steel at f=0.06 

Fig. 3.22 Variation of Roll Torque with R/h^ 72 

at different Reductions 
for Steel at f=0.06 

Fig. 3.23 Normal Pressure Distribution Along the 

Arc of Contact at Different Average 
coefficients of friction for Steel 
at R/h^=^^ ; %r= 16 

Fig. 3.24 Shear Stress Distribution Along the 

Arc of Contact at Different coefficients 
of Friction for Steel at R/h^=65 ; %r = ^6 

Fig. 3.25 Variation of Roll Separating Force with 75 ’ 

Coefficient of Friction for Steel 
at R/h^=^S ; %r= 

Fig. 3.26 Variation of Roll Torque with Coefficient 76 

xi 



Fig. 3.27 


Fie- 3.28 


Fig. 3.29 
Fig. 3.30 
Fig. 3.31 
Fig. 3.32 


Normal Preaaura diatribution Along the 
Arc of Contact Uithout and Uith Strain 
Hardening, (a) Minimum and (b) Maximum 
expected in the range 
Shear Streaa Diatribution Along the 
Arc of Contact Uithout and Uith Strain 
Hardening (a) Minimum and (b) Maximum 
expected in the range 
Normal Strain Rate Contoura 
Equivalent Strain Rate Contoura 
Deviatoric Streaa Invariant Contoura 
Equivalent Strain Contoura 



LIST OF NOTATIONS 


Synbol Deacrlption 

A Area of the Domain 

A Area of an Element 

B Matrix of Differential Coefficients of Shape 

Functions 

f Average Coefficient of Friction 

f|^ Boundary Elemental Matrix 

F Global Load Matrix 

h^ Initial Thickness of the Uork Piece 

Final Thickness of the Uork Piece 

H Matrix Containing Defferential Coefficients of 

Shape Functions, ( A Part of the Matrix B,used for 
strains ) 

jJ| Determinant of the Jacobian Matrix 

k^ Relative Sign of Boundary Tractions 

k Elemental Area Matrix 

kj^ Part of the Elemental Area Matrix Resulting from 

Momentum Equation 

k^ Part of the Elemental Area Matrix Resulting from 


xiii 



Penalization 


Global Stiffness Matrix 


Eleroental Coefficient Matrix for Strains 


Global Coefficient Matrix for Strains 


Length of the Inlet Zone 


Length of the Exit Zone 


in 

n 

nbelesi 

nelein 

n 

X 

n 


Multiplication Matrix 

Multiplication Factor for the Calculation of inlet 
and exit zone lengths 
Number of Boundary Elements 
Number of Domain Elements 
Unit Normal Component in X-Direction 

Unit Normal Component in Y-Direction 


n 


N^ 


N 


r 

R 

R 

S 


Unit Normal Vector 

One-Dimensional matrix of Biquadratic 

Shape Functions 

Matrix of Biquadratic Shape Functions for the 
Approximation of Velocities 
Hydrostatic Pressure 

Elemental Right Hand Side Vector for Strains 
Global Right Hand Side 
Roll Radius 

Deviatoric Stress Tensor 


xiv 



Deviatorlc Stress Components 

Deviatorlc Stress Invarient 

natrix of Boundary Elemental Nodal Tractions 
Normal Traction on the Uork Piece-Roll Interface 

Tangential Traction on the Uork Piece-Roll 

Interface 

Traction on the Uork Piece-Roll Interface in the 
X-Direct ion 

Traction on the Uork Piece-Roll Interface in the 
Y-Direction 

Velocity Component in X-Direction 
Resultant Velocity Normal to the Interface 

Known Velocity in X-Direction at the Inlet 

Known Velocity in X-Dlrection at the Exit 

Velocity Component in Y-Direction 

Known Velocity in Y-Direction at the Inlet 

Known Velocity in Y-Dlrectlon at the Exit 

Velocity Vector 

Global Matrix of Nodal Velocities 

Elemental Matrix of Nodal Velocities 
Global Matrix of Nodal Ueight Functions 

Matrix of ElementalNodal Ueight Functions for 


Strains 



Global Matrix of Nodal Height Functions for Strains 


X-Coordinate 

Y-Coordinate 

Angle Made By the Normal to the Interface 
Y-Axis at any point on the Interface 
Angle at a Global Node on the Interface^ k 

Neutral Angle 

Partial Derivative Operator 

The Vector Differential Operator , Del 

Strain Rate Tensor 

Strain Rate Components 

Deviatoric Part of the Strain Rate 

Equivalent Strain Rate 

Equivalent Strain 

The Boundary of the Domain 

Elemental Boundary 
The Penalty Factor 

Constant Used to Evaluate the Penalty Factor 
Density 

The Stress Tensor 


with 


The Stress Components 



Initial Yield Stress of the Metal 

Variable Yield Stress of the Metal 

Natural Coordinates for Domain Elements 
Proportionality Constant 

Constant Used for the Non-Dimensionalization of /j 
Natural Coordinates for Boundary Elements 

Boundary 

Corresponding to the Velocities 
Corresponding to the Strains or Tangential 


0 

n 

1 
2 
X 

y 

Superscripts 

be 

e 


xvii 


Initial Value 
Normal Component 
Inlet 
Exit 

In the Direction of X-Axis 
In the Direction of Y-Axis 

Boundary Element 
Domain Element 

Non-dimensionalized Quantity 


o 

ty 

y 

2tJ 



b 


V 

s 



CHAPTER 1 


INTRODUCTION 

1.1. General Introduction 

Rolling ia one of the moat important induatrial metalworking 

proceaaea becauae of ita extenaive uae in the production of a 

large variety of induatrial ateel products. Moreover, aa ia 

evident from the extent of published literature since the 

beginning of this century, rolling is apparently one of the moat 

investigated of all the known metal forming processes. 

1<2. Mechanics of Rolling Process 

In reality as the material approaches the rolls it first 

undergoes elastic deformation till it reaches the yield point and 

a large plastic deformation zone follows. As the material exits 

the deformation zone there will bean elastic spring back. However, 

the magnitude of the elastic strains is negligibly small compared 

to the large amounts of plastic deformation that the material 

undergoes during the rolling process. Typical orders of elastic 

-4 

and plastic strains usually found are 10 and 0.25 respectively. 
Therefore ignoring the elastic strains the material may be assumed 


to be 

rigid 

till it 

comes first in 

contact 

with the 

rolls 

and 

after 

it leaves the 

contact of the 

rolls 

Further 

during 

the 

deformation 

proceas 

in the deformation 

zone the 

strains 

are 


totally plastic. 

Many experimental and theoretical investigations have shown 
that the velocities, strains and stresses are not uniformly 
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distributed serosa the sections. The velocity across the sections 
of the work in the unstrained zone remains uniform while in the 
deformation zone the outer layers of the metal towards the rolls 
move faster than the inner layers. Near the entry the outer layers 
have a less velocity than the roll periphery. But the velocity of 
the outer layers (represented by V ) of the work piece increases 

£t 

as it moves towards the exit and eventually it becomes more than 
the roll peripheral velocity(represent ed by V^) as it reaches the 
exit. Although the transition occurs over a zone, because of the 
lack of exact data many researchers considered this as a point. 
The zone in which V is less than V is called forward slip region 
and in this region the friction force tends to retard the motion 
of the work piece. The point at which V^is equal to is called 
the neutral point. 

Analytical expressions for the determination of the neutral 
point have been developed by some researchers. Notable among them 
is the Avitzur’s formula [ 1 ]: 


ot 


n 



tan 



1 


- 1 )- 



h 


1 


) 



- 1 ) 


+ ).] 

(2/y3)o-^ 

The strain rates developed in the rolling process where the 
roll velocities are not very high are appreciably small. The 
entire work done by the tool in the plastic deformation of the 
workpiece is converted into heat leading to an increase in the 
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temperature of the workpiece. Strain-rates, temperatures and the 
strains that the material undergoes can affect its yield strength. 
However, the effect of the former two viz. strain-rate and 
temperature effects can be ignored in a cold rolling process since 
they are not substantially high whereas the strain hardening 
effect is substantial to be considered. 

Unlike most of the other metal forming processes where the 
tool remains stationary, in the flat rolling process the rolls 
will also be in movement thus the mathematical modelling of the 
process is difficult. The friction coefficient constantly changes 
from point to point along the arc of contact becoming zero at the 
neutral point. Ilany analytical expressions have been proposed in 
the literature to model the friction behaviour at the interface in 
the rolling process but none has predicted its variation 
correctly. 

For a successful design of rolling equipment and the process 
proper the important parameters to be considered are the normal 
pressure distribution along the interface, roll separating force 
and the roll torque. Analytical estimation of these parameters 
entails the analysis of the strainrates, strains in the workpiece 
during the deformation process. 

1.3. Previous Uork 

Many researchers carried out a variety of analytical and 
experimental investigations on the complex problem of rolling 
process . 

The pioneering work on the theory of rolling was done by Von 
Karman [ 2] in 1925, and then by Orowon [ 3] in 1943. As per the 
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review given in [ 5 ], von Karstan introduced for the first time 
the slab method of analysis where the assumption of planes remain 
was made, while the later work came out to be one of the most 
comprehensive of all analytical studies where Orowon accounted for 
all the various phenomena of rolling process and developed a 
differential equation. Orowon could not give an analytical 
solution to it but he presented a graphical solution. Considering 
the complexity of Orowon’s equations , the later researchers, 
especially Bland & Ford [ 4 } introduced some simplifying 

assumptions and developed analytical expressions. Alexandar[ 5 ] 

presented a comprehensive solution to Von Karman's original 
differential equation without taking into consideration the 
effects of inhomogeneity of deformation. Lahoti et al [ 6 ] 

combined the best of the theories then existing and solved for 
temperatures , roll separating force and roll torque with the aid of 
a computer-aided simulation and compared his results with the 
experiments in the literature. Numerous slip-line field and 
upper-bound solutions have also been proposed for the plain-strain 
rolling case. A few popular ones are the Alexander’s [ 7 ] 

slip-line field solution for hot rolling, the slip-line field 
solution of Firbank etal [ g ]for cold rolling and the only known 
upper bound solution of Jhonson etal [ 9 ]. The slip-line and 

upper bound techniques can give good solutions only to a limited 
number of cases and the material is invariably assumed to be 
perfectly plastic. The solution is also not unique in these cases. 

The advent of the application of the powerful technique of 
finite element method to metal forming processes is of recent ^ i 
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past and It Indeed paved a way to a great multitude of analytical 
solutions which till date have taken different turns towards the 
advancement in unraveling the problems involved in the rolling 
process, especially the neutral point phenomenon. The initial uses 
of FEH in metal forming problems were made by Ziekiewicz etal 
[ 10 ] and Dawson[ 11 ] who treated the rolling problem in terms 

of viscoplasticity where the metal is considered to flow as an 
incompressible non-Newtonian fluid with prescribed boundary 
velocities. Dawson also included the thermal analysis . G. J . Li . et 
al[ 12 ] have analyzed the problem of cold flat rolling by 

rigid-plastic finite element analysis by considering the friction 
relation developed by Chen and KobayashiC 13 ] . At the same time 
nori etal [14] analyzed in the same direction but incorporated 
slight compressibility of the work material. Yhu et al [ 15 ] 

accounted for the variable coefficient of friction by using the 
experimental data in the literature and obtained a better 
correlation with the experimental results. Richelson [ 16 ] 

studied the effect of various friction models on the distribution 
of the normal and tangential tractions in cold plane-strain 
rolling using rigid-plastic FEH. However, the final important 
parameters of roll force and roll torque were not estimated in 
this work rendering it to be incomplete. An interesting case of 
asymmetric cold rolling was dealt with by Pietrzyk [ 17 ] where 

the problem of plane-strain rolling was analyzed with different 
angular velocities and different diameters of the rolls 
considered. The friction coefficient was used differently for 

the two rolls in the pair.Hori etal [18] tried to incorporate 


5 



th« elastic deformation of the rolla by boundary element method 
while the strip was analyzed by rigid-plastic finite element 
method. It was illustrated that the normal pressure is lower with 
deformed rolla, than that with rigid rolls , in contrast to what 

was believed by many researchers. Some 3-D analyses have been 

conducted by researchers notable among them are [ 19 - 24 ] . 

A good number of experimental investigations have also been 
conducted in rolling. Al-Salehi et al [ 25 ] has presented a 

comprehensive set of experimental results by measuring the roll 

pressure distribution in the roll gap during cold metal rolling 

under plane-strain conditions for aluminium, copper and mild 
steel. The difference in pressure distribution with different R/hl 
values was highlighted. It was shown that neutral point occurs in 
all cases of rolling very close to the exit. Roll force and roll 
torque values were measured and presented for aluminium and copper 
at various reduction ratios. Shida etal f26J have measured the 
roll force and roll torque in case of cold rolling of mild steel. 
1.4. Objectives of the present thesis 

From the above study of the finite element literature it is 
clear that while the initial works ignored the strain-hardening 
effect, the later investigations except for [ 15 ], considered a 
frictional relation in which the so called friction factor ’m’ is 
used which does not seem to have any rational basis. Also few 
works have estimated the error which results from the assumption 
of perfectly plastic material. Some researchers did not determine 
the practically useful quantities like the roll force and roll 
torque while many others did not compare their results with any 
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experimental work. 

In the present work the important design parameters o£ the 
rolling process viz. roll force, roll torque and normal pressure^ 
shear stress distributions are estimated for various cases 
involving three metals under the plane-strain assumption and by 
considering the effect of strain hardening. For treating friction 
at the work piece-roll interface and to model the neutral point 
phenomenon, the traditional Coulomb's friction relation is used in 
which the average coefficient friction is taken from the same 
experimental works in the literature from which the analyzed cases 
are considered. The analytical results for roll force, roll torque 
and normal pressure distribution are compared with the 

experimental results. A detailed parametric study is carried out 
and the percentage error that is likely to accumulate with the 
ignorance of strain hardening effect is highlighted. 

Determination of the above mentioned design parameters of the 
rolling process needs the knowledge of stress and strain 
distributions in the work piece. These are estimated by treating 
the rolling problem as a plane-strain steady state problem using 
Eulerian formulation with strain rate as the measure of 
deformation. The material is assumed to be rigid-plastic strain 
hardening. Temperature and strain-rate effects are ignored as 
these are usually small. For simplicity, the rolls are considered 
to be rigid. The domain includes sufficient portions of the inlet 
and exit zones to take the advantage of uniform conditions 
prevailing there. Galerkin weighed residual technique coupled with 
the direct penalty method is used for weak formulation. The domain 
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Is discretized into snsll elements and biquadratic approximation 
for velocity components is used to convert the continuity and 
momentum equations into a set of nonlinear algebraic equations. 
These equations are solved by iteration. First a converged solution 
is obtained without the strain hardening effect. This solution is 
then used as the initial approximation to solve the system of 
equations for strain hardening problem. The knowledge of the 
equivalent strain required for updating the material properties at 
each iteration is obtained by integrating the equivalent strain 
rate. The converged solution for nodal velocities is used to 
estimate secondary quantities viz. strain rat es , stresses , rol 1 
separating force, roll torque and normal pressure and the shear 
stress distributions. 

Our modelling is such that the roll velocity la not required 
in our non-dimensional formulation. 

1.5. Scope of the Present Uork 

In this analysis, the material is assumed to be rigid plastic 
strain hardening which yields according to the von nises 
criterion. Although this assumption is valid in the deformation 
zone the elastic effects are significant at the entry and exit and 
therefore, should be included in the analysis. However, for 
simplicity, these effects are Ignored in present analysis. The 
rolls are assumed to be rigid. Roll flattening has a propensity to 
decrease the roll forces since the effective reduction with 
elastically deformed rolls is lesser [ 18 ]. 

The effect of temperature and strain rate effects on the 
yield strength of the material are ignored in this work because of 
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the two reasons, besides the fact that their inclusion renders the 
analysis to be too complex. The first reason is that it is 
strongly expected that whatever the direct effect of strain ^»rate 
on the yield strength is present it is marginalized by the 
temperature softening. Secondly the roll speeds are not so high as 
to increase the strain rates to magnitudes. It is observable that 
no previous work has considered both strain and strain rate 
effects simultaneously. 

The most important of all assumptions is that of average 
coefficient of friction. But in practice this is not the case 
where the coefficient of friction continuously varies along the 
interface, becoming zero at the neutral point. 

The roll velocity can enter the problem only through the 
dependence of the yield stress on the strain rate or through the 
roll velocity dependent friction relation. But since we are not 
incorporating any of these two effects, there is no dependence of 
various parameters on the roll velocity in this formulation. 

The present analysis investigates into the applicability of 
the traditional Coloumb's friction relation for the treatment of 
the friction at the workpiece-roll interface in the problem of 
cold flat plane-strain rolling of plates where the width of the 
plate is considerably large compared to its thickness, at steady 
state. In this relation the frictional stress is a constant 
multiple of the normal pressure, this constant being the average 
coefficient of friction and in this analysis this has been 
obtained from Al-Salehi’s[ 25 ] experimental results for aluminium 
and copper and from Shida etal (261 for steel. The deformation 
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of the plate la analyzed by the rigid-plaat ic finite element 
technique coupled with the direct penalty method and Eulerian 
formulation. An endeavour ia made in the analyaia to illuatrate 
the effect of atraln hardening on the important parametera of roll 
force, roll torque and the normal preaaure diatribution which are 
of utmoat importance in the mill deeign. The reaaonably good 
correlation between experimental reaulta [ 25,26 ] and theoretical 
reaulta over a range of input parametera ahowa that thia model haa 
certain validity and if extended with variable coefficient 
friction, non-rigid rolla and elaatic entry and exit it can 
produce more accurate reaulta. 

1.6. Plan of the Theaia 

The theaia ia organized aa followa. 

In the aecond chapter the mathematical modelling and the 
Galerkin finite element formulation of the problem of plane-atrain 
cold flat rolling are preaented. The application of the boundary 
equationa and the aolution procedure are alao diacuaaed in thia 
chapter . 

In chapter 3 the reaulta obtained by running the model for a 
range of caaea are preaented along with the compariaon with the 
experimental reaulta and a diacuaaion on parametric atudy. A 
concluaion ia drawn and auggeationa for further work in the field 
are given. 
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CHAPTER 2 


MATHEMATICAL MODELLING OF COLD PLANE-STRAIN ROLLING PROCESS 

2.1. Introduction 

In this chapter the matheinatical model for flat cold rolling of 
plates is developed. Considering that the width of the plate is 
more than ten times the initial thickness , the process is modelled 
as a steady state plane-strain .problem. The constitutive equations 
for the mechanical behaviour of the work material are stated. The 
interaction between the work and its surroundings are represented 
by appropriate boundary conditions. Non-dimensionalization of 
governing relations and Galerkin formulation are presented at the 
end of this chapter. 

2.2. Governing equations 

The mechanical behavior of the metal is governed by the 
following two equations. 

(1) Conservation of mass or continuity equation : 

p + p div.y =0 (2.1) 

(2) Conservation of momentum or momentum equation : 

p V = (for zero body forces) (2.2) 

where the dot above V denotes the material time derivative. 
The constitutive relations , besides these governing equations , that 
are used to model the metal behavior are as follows. 

In flow formulation the measure of deformation is the rate 
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of strain tensor, which stay mathemat ical ly be expressed as 


IT 1 ? 

or ‘ij= i Ciij* 5ip C2.3 

In order to express the stress as a function of strain rate 
in a convenient form the stress and strain rate tensors can be 
divided into two parta_ 


f= -P 1+ f or (2.4) 

where p=-|tr <y is the hydrostatic part, S is the deviatoric part 
and is the kronecker delta. Similarly, 


= |tr 


•f jD 


or 


ij 


'rkk^^ij 


+ £ 


ij 


(2.5) 


where £ is the deviatoric part of the strain— * rate tensor. In 
plastic deformation, since there is no change in volume the 
hydrostatic part is not related to the deformation. Another 
consequence of volume constancy is that the hydrostatic part of 
the strain rate is zero and its deviatoric part is the strain rate 
tensor itself. Then the deviatoric parts of stress and strain rate 
tensors can be related as : 


5 “ 2^^ or S j * ^fjt£ ^ j ( 2 . f ) 

where 2(j is the proportionality constant. Further if we define 
the second invariant of deviatoric stress as 




and second invariant of strain rate 


(also known as 


(2.7) 

equivalent 
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atr&ln rate ) aa 





then it followa that 


( 2 . 8 ) 


S=iJ£ . (2.9) 

For inetala which yield according to von Miaea criterion. 


Cf 

S= 

^ ya 


( 2 . 10 ) 


and p for such inetala comes out to be 

"y i C2.11) 

where a is the yield strength of the metal. In general cf may be 

y* y 

a function of equivalent strain , temperature and equivalent strain 
rate. However, if the deformation process is such that the 
temperature and viaco-plaatic effects on the mechanical properties 
can be neglected then it is a function of only equivalent strain. 

For a strain hardening metal the yield stress is related to 
equivalent strain by an equation. 


o =o- (l+Ki)*' (2.12) 

y o 

where ^ is yield stress of the metal at zero plastic strain and K 
and n are metal dependent coefficients determined from 
experiments. The equivalent strain is obtained by time integration 
of equivalent strain rate as given below! 

t 

2=J J dt (2.13) 

0 

2.3. Equations for steady-state two-dimensional problem 
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The work 

is assumed to 

be 

very wide 

as compared 

to 

its 

thickness and 

consequently 

the 

effect of 

the spread 

may 

be 


neglected and the problem may be considered as plane-strain. The 
rolls are assumed to be rigid and the metal is assumed to be 
rigid-plastic strain hardening yielding according to von Mises 
criterion. The temperature and visco-plastic effects are neglected 
because it has been observed that the strain rates in the cases of 
cold rolling process considered are ignorably small. Only half of 
the domain is considered for analysis because of symmetry, see 

figure 2.1, and the rectangular cartesian coordinate system as 

1 

depicted in figure 2.2 is used. 

For plane-strain problem the velocity vector has only two 
components , 

V=ui+vj (2.14) 

while the stress and strain-rate tensors have the following form: 



o* 

XX 

o* 

xy 


lO 

XX 

xy 


<y 

L yx 

O' 

yy 


yx 

i) 

yy J 


t 


The process is analyzed at steady state and therefore the 
governing equations as reduced to 2-D are: 


du 

Sx 



(2.16) 


/■ du^ du -v 

c 3 


^ ss ss 

^2 + ** + — 

dx dx dy 


r SV . dV .. 

C 3 


* . 3S 

+ __y* + __yy 

Sy Sx Sy 


(2.17) 


(2.18) 


where 
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S =2p" 

XX XX ^x 


S =2^£ =2pt- 

yy yy ^y 


S =2/J£ 
xy xy 


C ) 


(2.19) 



2.4. Boundary conditions 

The boundary conditions specified on any boundary can be 
classified as essential boundary conditions i.e., the 
specification of the velocity vector V, natural boundary 
conditions i.e., the specification of the traction vector t~<y.n or 
mixed boundary conditions i . e. .specif ication of one velocity 
component and one traction component. For the problem under 
consideration the boundary conditions are as f ollovs ( f igure 2.3). 

(1) Entry and exit boundaries (AFaDE): 

The control volume is so selected that its entry and exit 
boundaries are sufficiently far away from the either side of the 
deformation zone. This ensures a uniform variation of the velocity 
before and after deformation and allows the following 
boundarycondit ions for consideration on these regions! 


II 

and 

o 

11 

> 

on side 

DE 

u=Uj 

and 

o 

II 

> 

. on side 

AF 


(2.23) 


i.e., all the material particles on these boundaries have the 
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preacribed velocitiea and at entry and exit reapectively , 
and the velocitiea of theae partlclea in Y-dlrection are zero. 

(2) The top free aurfacea (AB&CD): 

Since theae aurfacea are traction free, 

t =0 and t =0 
X y 

Alao aince the flow ia only in the X-direct ion , v=0 . But only 
one condition out of the two (t =0,v=0) haa to be choaen. So the 

y 

final choice of the boundary conditiona ia 

t^=0 and v=0 (2.24) 

(3) The axia of ayfflinetry(FE) : 

Becauae of aymmetry, v and t^ are zero on thia region of the 
boundary. Thua 

v=0 and t^=0 on the axia of ayminetry. (2.25) 

(4) Uork-roll interface (B©): 

The component of velocity in the direction normal to the 

roll-plate interface at any point on the interface la zero. 

u =0 or u+v tana=0 (2.26) 

n 

where u and v are the velocities in the directions x and y at any 
point on the interface and a. is the angle made by the unit normal 
to the Interface at that point with the y-axia. At the first 
contact point on the interface towards entry the angle ot is equal 
to the angle of bite. 

The complex and Important part of the analysis of rolling 
process la the treatment of the work-roll interface friction and 
modelling of the neutral point phenomenon. Many researchers have 
tried to tackle this problem by using rather arbitrary functional 
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relations for shear stress in such a way that the position of 
neutral point is constrained to be within the inlet and exit of 
the interface, in accordance with the experimental findings. 
According to Al-Salehi etal,[ 25 ] the neutral point, such that 
the coefficient of friction falls to zero and a shear stress 
reversal occurs, is present in every case and close to the exit 
plane. However, none of the works ha^ tried to incorporate the 
experimentally measured coefficient of friction either directly or 
Indirectly. Such an endeavor is made in this analysis and thus the 
following relation is used: 

I „l (2 

where f is the coefficient of friction. In this analysis an 

% * 

average value of fis used assuming that it remains constant 
throughout the deformation process. These values for individual 
cases considered are taken as reported by Al-Salehl etal . 

Further the shear stress in the eq. (2.27) Is subject to the 
constraint 

</ cr 

if t > then t = -^ 

® yi ® ya 

The method of application of this boundary condition is explained 
later in the next chapter in detail. 

2.5. Non-dimensionalization 

Non-dimenslonallzat ion is usually done to avoid 
ill-conditioning and numerical difficulties. The 

non-dimensionalization of all physical quantities using 
characteristic dimensional quantities is presented in the 
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following lines: 


rw 

x= 


X 




u=- 


u 


U 


ft 


(2.28) 

*^ft 


I f we use 


o 

^o y3(2U~7h") ' 


(2.29) 


to non-dimenaionalize (j then 


the expression for may- 


be 


written 


as , 



(2 . 30) 


where the non-dimensionalized strain rate invariant is given by 



(2.31) 


For a perfectly plastic material 


P = lA 


(2.32) 


The non-dimensionalization of the other quantities is as 
follows : 


s s 

xy xy 

^xy="^;2fl2/h""‘<5;7y3" 


XX 


XX 


“XX 

s =— _-yy— =— -Yy— 

yy ^02^2^2 


(2.33) 


The non-dimenslonalizatlon of the continuity and momentum 


IS 



equatlona la obtained by aubatltutlng the eq8.(2.28) to (2.33) in 
eqa.C2.16) to (2.18) . 


dx Sy 


(2.34) 


pu h ^ ~ ~ ^ ~ ^ ss as - 

_f __¥* + __XY Vo 

^ ^x Sy ' Sx ^ dx Sy ' 


alinilarly , 

pU.,h 


r~ av ~ 

|u --+V 

av 

V -[ 

.as .. 

__Y? + __Yy 1 

^ ax 

ay 

i ay ' 

■ ax ay ^ 


]-° 


(2.35) 


(2.36) 


Since the Reynolda number 



ia found to be typically of 


— g 

the order of 10 for all the caaea conaidered, the inertial terms 
may be ignored. Thua the summarised non-dimensional forms of the 
continuity and momentum equations are: 





^x ^y 


(2.37) 


- ^2 ^ 


Sx 


, as 

f __xx^ xy 

^ ax ay 



(2.38) 


_ fe + r __Y5+__YY 1=0 (2.39) 

Jy ^ 49 x ay * 


2.6. Galerkin Finite Element Formulation 

2 . 6 . 1 .Galerkin Ueak Formulation with Direct Penalty Method: 
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In thla analyaia the Galerkin weighed reaidual method 
coupled with direct penalization ia uaed for the finite element 
formulation of the governing equationa. In the mixed 
or preaaure-velocity formulation where preaaure ia alao a nodal 
variable, the coefficient matrix becomea a non-poaitive-def inite 
matrix becauae of the zeroea appearing along the principal 
diagonal, therefore entailing a aolution method employing complete 
pivoting. Alao aince the preaaure doea not occur aa a degree of 
freedom at every node, the computer implementation of aaaembly 
becomea more complicated [ 27 ]. To overcome theae problema and 
becauae of aome of the other advantagea with direct penalty 
formulation with reduced integration [ 2S ], thia method ia 
adopted in the preaent analyaia. 

In direct penalty formulation the preaaure in the momentum 
equationa ia eliminated uaing the following relationahip 


where 



(2.40) 


and X, which ia called the penalty factor, ia a very large 
number. In the limit aa X -+ oo , equation (2.40) reducea to equation 
(2. 37). The value of X ia uaually choaen in auch a way that 


the upper limit being dependent on the computer. 

The weak form ia derived in thia aub-aection and the finite 
element formulation ia preaented in the next aub-aection. Let u 
and V be the functiona which aatiafy all the eaaential boundary 
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condltione exactly. Then u and v constitute a weak solution if the 
following integral equation is satisfied. 


J 


Sx ^ ay ^ 



ss 

ss 


ss 



+ __?Y 

^-v^ 

- 11 

dA =0 

Sx 

Sy 

^ ^ Sx 

Sy 



(2.41) 


where w^ and w^ are weight functions which satisfy the homogeneous 
versions of the boundary conditions and 'A’ represents the area of 
the domain. 

Performing the above integration by parts the following weak 
form can be obtained: 


I 


Sw 

S\j 

, dw 


r Sv Sv 1 


1 

1 

1- 

1 

! 

\< 

mm* 

hf s 

4. Q 

I 

I 

m 

h 

I 

I 

4. Q 

Sx 

Sy J 

1 XX 

Sx 

xy 

~ ^ ~ 

L dy Sx j 

T ^ 

yy ^y 


dA 


I I fw S n + w S n +w S n +w S n "I - fw pn +w pn "1 I dr =0 
J u XX X u xy y v yx x v yy y-^ v u x y-'J 

r (2.42) 

where T represents the boundary of the entire domain under 

\ 

consideration . 

Substituting for p from eq. (2.40) and for form 

eq. (2.19) only in the area integral, the relation in eq . (2.4) 

only in the boundary integral and using the relation 


'' y-'L yx yy-*^ y-' 


(2 . 43) 
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we obtain the final weak form ae 


Sv dw 


J Sx ay ^ ^ Sx Sv ^ '- 


Sx Sy 


Sv ~ 

|4C ---+jC 

XX xy 



- ft w +t w 1 dr =0 
J '' X u y 


(2.44) 


2.6.2 Finite Element Approximation: 

In contrast to the traditional Galerkin-Ueighed Residual 
Method where a global approximation is chosen for the field 
variable, in finite element method the domain is discretized into 
small elements and the above integral is evaluated over each 
element after suitably selecting the approximation for the field 
variable over each element. Then the elemental expressions are 
assembled to yield the global FEM expression. 

The eq. (2.44) suggests that the convergence criteria will be 
satisfied if u and v are chosen to be bilinear functions of x and 
y. So one can use bilinear interpolation for u and v. 
However , higher polynomials give better results. lie choose 
nine-noded Lagrangian element and corresponding biquadratic nodal 
shape functions for u and v, figure 2.4* 

The approximation for u and v is as follows: 
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( 2 . 45 ) 



V ]= [ NJ[ V*] 

In Galerkln forntulation the weight functions are approximated 
using the shape functions as that for velocity. 


I 

The vectors [v-l and [“*] respectively contain the velocities 
at the nodes of elements and their corresponding weight functions. 

The expressions for the shape functions for nine noded 
Lagrangian quadratic rectangular element are given by: 


“4=1 p'-O [»'*»] = 


“5=1 

Vl P'-OP-’*'] 




“ 6=1 


N,= 


where ? and rj are the local coordinates and are related to the 
global coordinates x and y with subparametr ic formulation as 
follows : 
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(2 . 48) 


^ 4 ^ 

x=E 


whert 






4 


y=E 


'‘2'iCi*0(i-«D 

C2.49) 


The boundary approximation, figure 2.5, is made consistent 

with the area approximation as follows. 



[ NJ [ “b*] 


(2.50) 


where [“bl is the vector of Lagrangian quadratic shape functions 
and is the vector of local boundary weight functions for a 

typical ' boundary element. Further we approximate t^ and t^ 

as , 



(2.51) 
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b 6 

where -{tj^ } ia the vector of local boundary tractlona for a 
typical boundary element. 

2.7. Finite Element Equations in Local Variables 

The following definitions of some vectors should precede the 
presentation of elemental expressions and their assembly, for 
compactness and ease of understanding. 

The strain-rate vector is defines as follows. 


or further 



C2.52) 


(2.53) 
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and 


where 


u 

~Zi 

Sx 

Sv 


ay 

^ ^Sv av . 

if— “♦— * 




= [ B ] [ U-| 


(2.54) 


[ B ] 


1 

y2 


Sx 

0 

ay 


0 

au^ 

ay 

1 '"'i 


an, 

4 

ax 

0 


1 "2 
-/2 ~~ 


0 

fN2 

ay 

1 ^^2 


■/2 


ax 


au^ 

ax 

0 


an 

1 __9 


0 

dNj 

ay 

t 

dx 


(2.55) 

from 


Substituting the the area elemental approximation 
eq . (2.45) to eq. (2.49) in the area terms of eq. (2.44) and 
noting that 


du^^v 

ax ay 


= [«]^[B][V'] 


and 


av av 
ax ay 


where 


E »] 


1 

1 

0 


.we obtain the elemental area term in matrix form as 
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A* 

|22[ U*]^[ B ]^[ B ][ v']|jdA* (2.56) 

Further aubstltutlne the boundary approximation from eqs . (2.50) 

and (2.51) in the boundary term of eq . (2.44) we obtain the 

elemental boundary term as 


r^e 


(2.57) 


To convert the above integral equations from the global 
coordinates to local coordinates , x and y are to be changed to ? 
and q using the following relation. 

+ 1 +1 

|( ) dx dy = j J ( )|J| d? dr> (2.58) 

-1 -1 


where |J|,the Jacobian, is the determinant of the Jacobian matrix, 


^x dx 

d? dy> 

Thus the eq. (2.56) and eq. (2.57) can be written in terms of 
natural coordinates as 
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+ 1 +1 


-1 -1 


and 


+ 1 +1 

J J {22[“*]^[B]^[B][n}]|Jl ■«: 

-1 -1 


(2.59) 


+ 1 

C2.60 

-1 

respectively where |jj is the Jacobian for the area terms and l^j^j 
is the Jacobian for the boundary term and C is the natural 
coordinate for the boundary linear elements. 

2.8 Global Assembly 

The final finite element equation is obtained by 
assembline the elemental area and boundary matrices into 
global area and boundary terms , respectively , and then equating as 
follows . 
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nel era 
E 

e=l 


+ 1 +1 

M ® *" 3C ® ]^C B ][ v*]J|j| de d77+ 

-1 -1 

+1 +1 

J I {2^[ B ]’'[ B ] [ V*]||J| K d„ 

-1 -1 

+ 1 

nbelemf r .. a« 

' 1 Jt C “b]^[ 1 


Further It can also be written aa 


nel era 
E 

e=l 




nbelem 

E 

be=l 




where 


+ 1 +1 

[ It 'll" I J {2?[ B ]^[ B ]}|J| d<f dfl 

-1 -1 


+ 1 +1 

[>'*]2= J J{[B]^[m][.]^[B]]|J| d« d„ 

-1 -1 

+ 1 

[ 'b“3= J{[ « 

-1 


bl 

C2.61) 


C2.62) 


(2.63) 


(2.64) 


(2.65) 


( 2 . 66 ) 


29 



The elemental area matrix is prepared by evaluating the 

integral by reduced integration technique where [k®] 3^ ia 
integrated by 3*3 gauss quadrature and [k®] 2 by 2*2 gauss 
quadrature. This ensures that j[k®J 2 ia singular so that even for 
fairly large values of X,the penalty factor, a non-trivial solution 
may be obtained. 

The assembly of the elemental area matrices into the global 
area matrix is done by transferring the elements corresponding to 
a local degree of freedom in each elemental area matrix to 
positions of corresponding global degrees of freedom in the global 
area matrix. Similar procedure is followed for the assembly of 
global boundary matrix. The assembled finite element equation 
may be written as 

[U]^[IC][V] = [U]^[F] C 2 .<S 7 ) 

where [ U ]| contains global nodal weight functions, [ ^ 3 

global stiffness matrix^ [ ^ 3 <^<^i^bains the global nodal velocities 

and [| F ]] is the global load matrix. 

Since [ U J is arbitrary we can eliminate it from both sides 
and thus we obtain the final general equation of FEH, 


[ K ][ V ] = [ F ] (2.68) 

2 . 9 . Application of Boundary conditions 

The application of natural and essential boundary conditions 
is done as follows. 

The natural boundary condition on the work-roll interface can 
be applied without requiring to evaluate the matrix [f^®] in the 
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following manner (figure 2.6). 

Let ua consider a typical boundary element on the work 

plece—roll Interface. The elemental boundary matrix for this 
element would be as given by eq. (2.66). At a middle node 'i' o£ 
tbia boundary element J, whose corresponding global node 

number la 'k' we can write the contribution of this node to the 
local boundary matrix and there by to the global load matrix as 
follows . 

+ 1 

C M2ic-r[ ‘b ]2i-r ft »'J^[ > 21-1 

>■ j 

-1 

+ 1 

tMlk “[‘bill = [ J[*j^[«b]lil < ^b> 21 

-1 

where [ F 32k-l ^ 1 2k elements of the global load 

matrix containing the odd and even net force contributions from 
the global boundary node 'k, respectively. Now for this node the 
relation in equation eq. (2.27) can be written as 
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Convtrtlne the and (ty)^^ Into aatri ce« analogous to, but 1 - 

boundary elemental traction matrix and multiplying both aides with 
+ 1 

T 

Jc ^ where [ N 2 N 3 ] ( 2 . 68 ) 

-1 

and replacing and (ty)j^ with { t^^ 2 ^_^and •lt^2i 

respectively, the following expression is obtained_ 

+ 1 

[J[nJ[NJ |JJ dC { t^} [coso.j^-fk^sinc*J 

-1 

+ 1 

-1 

This gives the contribution of one node of the element j'. 

Since a does not vary much from one local node to 

other in an element , similar evaluation for other two nodes and 
then summation over all three nodes of the element gives 

[ 2i-l ^b3 2i[®^*'“k'^^^8‘=°®“k3 ■ 

Further 

[ '■ 3.2k-l = [ f Izk 

Thus instead of evaluating [ F 32k-l t ^ 3 2k 

typical global Interface node k, we perform the following the matrix 
operations in order to apply this natural boundary condition. 

For each boundary node 'k' on the interface, 

(i) replace each element of (2k-l)th row of global stiffness 
matrix by: t 


32 



[ element of (2k-l)th row « (cosa -fk eina )- 

K fl k 

element of (2k)th row in the same column « (sina, +fk coaa. )1 

k a k -I 

(11) replace (2k~l)th element in global load matrix by zero. 

The boundary condition expreaaed in eq. (2.26) at a typical 
global node 'i' la applied through the following atepa. 

(1) Replace all the elementa of (2k)th row of the global 
etiffneaa matrix and that in the global load matrix by zeroea. 

(li) Subatitute '1’ in the elemental poaition of C2k)th row and 
(2k-l)th column, and the value of tana, evaluated with a at node i^ 
in the elemental poaition of (2k)th row and (2k)th column of the 
global atiffneaa matrix. 

2.10. Formulation for Strain Hardening 

The equivalent atrain may be obtained by the time integration 
of equivalent atrain ratea aa , 




£ dt 


(2.72) 


where £ and £ are the equivalent atrain and atrain ratea 
reapectively and the factor y| ia introduced ao that the 
equivalent atrain ia equivalent to the axial atrain in a tenaile 
teat. Differentiation of the above equation leada to 


d£ 1 - 
dt ‘ya ^ 


(2.73) 


where ^ria the material time derivative aa it ia in Eulerian 
dt 

formulation. For 2-D ateady atate ca3e,aB oura, the above equation 
can be written in non-dimenaionallzed form aa 







(2.74) 


33 



This equation has to be solved along with the condition ^=0 at the 
inlet . 

Galerkin foroiulation of the above equation over an element 
gives the following integral: 


A* 


C2.75) 


where w^ ia the weight function. 

The integral may be expreaaed in matrix form aa given below. 


I 


[ u V] 


Sx 

S£ 

Ai/ 


- w 


8 ya 


3 ^ 


dx dy=0 


(2.76) 


In the finite element approximation for equivalent atraina and 
the weight function the aame ahape functiona aa uaed in the 
velocity approximation are uaed. 


2=£ N 2 = i 
1 




NlH 


w^ 


(2.77) 


where <{ N f is the 1~D array of quadratic shape functions as 
defined below, 

i N‘f = [Nj N 2 ^9 3 (2.78) 

and the differentiation of atraina in equation (2.77) leada to 
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(2.79) 


Sc 

Sx 

sZ 

Sy 


[ H 2^ 


where [ H ] 


is defined &s given below. 




!“2 


dN 

— 

0 


0 ... 

— 

Sx 


dx 


dx 

0 


0 

f”2 

0 




dy 



0 

su^ 

Sy 

(2.80) 


Substituting the approximation for equivalent strain , weight 
function and velocities in eqs . (2.77) and (2.45) we obtain the 
following integral equation. 


nelem 

s: 

e=l 



k1]^[ V*]’^[ NJ^[ h 


dx dy- 


y||( [ 


= 0 


(2.81) 


To facilitate the numerical integration, the variables are changed 
to the natural coordinates. The final expression can be obtained 
as 



nelem 


£ 

e=l 


([ 


nel em 


[' ■I) 


(2.82) 
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where 


+ 1 +1 

[ O ' I J {[ NJ^[ H ]} |J| dlf d„ C2.83) 

-1 -1 

ia the elemental coefficient matrix for the evaluation of 
equivalent atraina and 

+ 1 +1 

[ r®] = y| J J[[ N^]^J j |J| d-n (2.84) 

-1 -1 

T 

is the elemental right side vector. Since £ U^j is arbitrary, we 
obtain the following equation in terms of global matrices, 

[ ^ = [ ® ] (2.85) 

where the is the global coefficient matrix, ia the global 

nodal equivalent strains matrix and w is the global right hand 
side, for the evaluation of equivalent atraina. In equation (2.84) 

£ is the non-dimensional equivalent strain-rate which may directly 
be obtained at each gauss-point from the velocity field of 
previous Iteration. 

2.11. Implementation 

Uhile correct modelling of the physical problem into a 
mathematical equation forms the core of the entire analysis, 
selection of degree of refinement of mesh, value of the penalty 
factor, the convergence criteria, the threshold value of ju and 
other mesh related parameters decides the achievement of 
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reaaonable reaulte. It haa been experienced that the judgment 
about the later mentioned aelection of various parameters needs a 
vide knowledge of capabilities and liabilities of computers as 
well it entails a good deal of numerical experiments before a 
particular set of values la used for final results. A brief 
discussion la presented in the following lines about some of these 
factors . 

The length of the inlet and exit zones is usually taken as 
follows : 


1 and l^=nh 2/2 

where 1^ and 1^ are the inlet and exit zone lengths respectively, 
'n' is a multiplying factor and inlet and exit 

thicknesses of the plate. The value for ‘n^ has to be selected such 
that uniform conditions prevail at the beginning of the inlet zone 
and end of the exit. A value of 6 has been selected in this 
analysis after conducting a few numerical experiments. 

The value of penalty factor X is usually selected as per the 
directions given elsewhere in this chapter. For this analysis X is 
found as follows: 

X=X/j 

*^0 

~ 2 5 

where X is a numerical value in the range 10 to 10 . Although the 

pressure remains fairly unaltered for a reasonably wide range of 

values of X, very high values may render the solution to be 

spurious. A value of X=10 has been chosen in this analysis and 

used constantly for all the cases. 
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I*'* this &n&lyai£i of rolling the neutr&l point is not known & 
priori. Therefore, for the first iteration an approximate neutral 
point is determined by the following Avitzur'a analytical formula: 



(27/3)^" I 


where 

a. = the neutral ancle 
n 

R = roll radius 

hi h 2 = the initial and final thicknesses of the plate 
m= the friction factor 

*^xb ^xf ” back and front tensions respectively 

O' = the initial yield stress of the metal 
Thus the relative signs k^ at each interface node is 
determined as follows: 


t 

8 

is 
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ve 

in 
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backward 

slip 
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k is (-) in the 
s ' 
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is 
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ve 
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k is (-) in the 
s 

n 

is 

(-) 

ve 

in 

the 

forward 

slip 
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forward slip zone 


This approximation is used only for the first iteration and the 
second iteration onwards the relative signs are determined from 
the computations of t^ and t^ from the solution of 
the immediately preceding iteration (the computation of t^ 
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converged solution 


and la explained later). 

The strategy followed for obtaining the 
with and without strain hardening is as follows. First a converged 
solution (nodal velocities) without strain hardening has been 
obtained. For this portion of the analysis the p is simply as 
given in the eq. (2.32), This solution is then used as the initial 
solution for the further analysis and the yield strength is 
updated by evaluating strains from immediately preceding solution 
in every successive iteration. Thus the value of ^ is evaluated 
from the eq . (2.31). 

Preliminary runs have been conducted with coarser mesh using 
various solution algorithms viz. gauss elimination method, IMSL 
routine with and without iterative solving, NAG routine etc. on the 
HP super mini computer to check the validity of our gauss 
elimination subroutine. Then final runs with refined mesh are 
conducted on CONVEX mini super computer to harness the high speed 
parallel computing using this solution routine. Because of the 
assumption of rigid portions in the inlet and exit zones, the high 
value of fj renders the coefficient matrix to be ill-conditioned, 
especially during strain hardening because the /j is updated. To 
counter this problem, ju is to be restricted to a threshold 
minimum, which is equal to the reciprocal of set off value of 
non-dimensional equivalent strain-rate. 

The convergence criteria used for velocity are as follows. 
All values of velocities less than 10~^ are considered to be equal 

to zero and are not checked for convergence. A percentage 

— /I ““2 

difference less than 10~ for u-velocitles and 10 for 
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v—vtloci 1 1 fifl b6tw66n fitoliit Ions of two successive iter&'tlons Is 
considered to mark the convergence of the solution , i.e., nodal 

velocities. This is done because the v-velocities are about two 
orders magnitude smaller than u-velocitles . 

2.12. Evaluation of the Secondary Quantities 
As explained above the solution to the problem is obtained in the 
form of nodal velocities. Thereby the secondary quantities viz. 
strain rates, deviatoric stresses, hydrostatic pressure and 
Interface boundary tractions are evaluated in each iteration to 
determine the non-linear parameters like ^,k etc., as well as at 

£t 

the end from the final converged solution for parametric study. 
The details of these calculations are explained in the following 
lines briefly. 

(I) The Hydrostatic Pressure at Gauss Points: 

The pressure at the gauss points is post-evaluated from the 
Lagrangian multiplier as follows: 

p=-xr ^^+^1 ] = - X 4 m B ]■{ 

^ dx dy 

This is evaluated at 2x2 gauss points. 

(II) Interface Boundary Tractions: 

Just like pressure the strain rates and the total stresses are 
also evaluated at 2x2 gauss points. The reciprocal of the 
non-dimensional strain-rate invariant helps in calculation of (j . 
The stresses at gauss points are linearly extrapolated on to the 
boundary nodes in a few steps. Then the tangential and normal 
tractions are evaluated at each interface node using the value of 
angle oi at that node. This is done in each iteration to find the 
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values of fj & as well as at the end for parametric study. 

(iil) The Roll Separating Force and Roll Torque: 

The roll separating force and torque are determined after 
obtaining a converged solution as follows. 

(a) Roll separating force: 

The roll separating force per unit width of the plate is 
evaluated as the integral of the vertical traction component (t ) 

y* 

over the arc of contact, expressed in terms of MN/m. Thus 

1 

s 

Roll separating force F =J t ds 

y y 

0 

where 1^ is the total length of the arc of contact. However, in 
this analysis, this integral la evaluated numerically using 2 
point gauss quadrature_ 

2 

‘y>i '•’I “i5 

where t is the normal traction at the gauss point *i,jJj is the 

y 

Jacobian of the transformation from the boundary coordinates to 
the natural coordinates and w^ is the weight function. 

(b) Roll torque: 

The roll torque is evaluated by numerically integrating, by 2 
point gauss quadrature , the elemental torque along the arc of 
contact as given below: 
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Roll torque =J Rt^ da =R 

0 

where " R " la the roll radiua,” t^” la the tangential traction 
at the gauaa point and |Jj and w. are aa defined above. 

(Iv) The Contoura: 

Contours are plotted to a how the dlatrlbutlon of certain 
aecondary parameters across the domain, namely for 2 ,£ and S. The 
contours for £ are plotted using the data at the nodes and for the 
later two, the plotting is done using the values evaluated at 2^2 
gauss points of the elements from the converged solution. 
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v = o 



Fie- 2.3 The Boundary of the Domain and the Boundary Conditions 



Element 
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Fifi. 2.4 A Typical Finite Element 
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Boundary Element tor Velocity 


Boundary Element tor 
Coordinate Trons-formation 


Fig. 2.5 A Typical Boundary Element 
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Fig. 2.6 The Boundary Tractions 




Fig. 2.7 The Finite Element Mesh 
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CHAPTER 3 


RESULTS AND DISCUSSION 

The finite element modelling of the physical problem of cold 
flat rolling developed in this work has been applied to a number 
of cases involving three metals and various non-dimensional input 
variables to illustrate its applicability. The metals considered 
are aluminium, mild steel and copper whose stress-strain relations 
are shown in table 3.1. The mat erials , dimensions and the 

strain-strain relations are considered in such a way that they 
are consistent with the cases considered in the works 
[12,15,25,26 ] so that the comparison is facilitated. In this work 
only comparisons with the experimental results are presented and 
rest of the emphasis is put on presenting a detailed parametric 
study and the effect of strain hardening. Comparison with the 
other analytical results is not given. Thus, in this chapter, the 
analytically obtained results are presented in four sections. In 
the first section the analytical results are compared with the 
experimental results of Al-Salehi et al [ 25 ] in case of 
aluminium and copper and Shida et al [ 26 ] in case of steel 
to validate the modelling and computations. In the second 
section, the parametric study, the effect of the three important 
rolling process variables viz. reduction ratio, the homogeneity 
factor (R/hl) and the coefficient of friction on the rolling 
parameters viz. normal pressure and shear stress distributions 
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alone the arc of contact, roll aeparating force and roll torque 
are discuaaed. In the third aectlon of thia chapter, the error 
that can enter by Ignorine the atrain hardening effect ia brought 
into light by comparing all the important deaign parametera 
without atrain hardening and with atrain hardening. The fourth 
aectlon marka the end of the preaentation of the reaulta where a 
atudy of contoura of normal at'Qiin ««ratea , equivalent atrain 
ratea , devlatorlc atreaaea and equivalent atrain ia preaented. In 
the penultimate aectlon of thia chapter concluaiona are drawn as 
a poat-analysla and evaluation of the work. In the last aectlon 
auggeationa for further research are included. 

3.1. Validation 

The validation of the reaulta of the variation of the roll 
force (per unit width) with reduction ratio is preaented in 
figures 3.1 to 3.6. The figures 3.7 to 3.12 show the comparison 
of the variation of the roll torque (per unit width) with 
reduction ratio. Finally the comparison of the distribution of 
pressure along the arc of contact is presented in figures 3.13 to 
3.16. 

It ia clear from the figa.( 3.1 ) to ( 3.6 ) that there is a 
good correlation between experimental and theoretical reaulta of 
the roll force in case of steel and aluminium and to some extent 
that for copper. However, at higher reductions the deviation 
between them increases except for the case of R/h^^ =12.5 for 
aluminium. Similarly as can be seen in figures 3.7 to 3.12, in the 
case of roll torque, except for the case of aluminium with 
R/h^=39, in all caaes the deviation between analytical and 
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experimental reaulta Increaaea. Uhlle it ia reported that the 
deviation ia becauae of the aaaumption of constant coefficient of 
friction [ 15 ],it can alao be attributed to the fact that at 
large valuea of reduction it becomea difficult to predict the 
neutral point accurately , which haa a direct effect on the roll 
torque. Although the effect may not be at par with that of 
friction, the rigid roll aaaumption and the ignorance of entry 
and exit elaatic effecta are alao expected to have contributed to 
the deviation of analytical and experimental reaulta. Figures 
3-13 to 3.16 show the comparison of analytical and experimental 
normal pressure diatribution along the arc of contact. Although 
there ia a considerable variation in magnitude, the important 
feature of multiple peaks of normal pressure at low values of 
R/hj^ and the friction-hill type of pressure diatribution at 
higher valuea have been reproduced. The multiple peaks at low 
R/hj^ have been predicted alao by the slip-line field solution 
of Firbank et al [ 8 ]. This suggests that the present model 
would be able to predict correct magnitudes of normal pressure 
over some range of valuea of R/h^ reasonably accurately. 

3.2. Parametric Study 

In this section a discussion on the effect of process 
variables viz. reduction ratio, R/h^^ and the coefficient of 
friction on the normal pressure distribution, shear stress 
distr ibut ion , rol 1 separating force and the roll torque la 
presented . 

3.2.1. Effect of reduction ratio and R/hj^ 

The figure 3.17 depicts the variation of normal preaaure<t^) 
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alone the arc of contact at different reductions. Uhlle the 
distribution is more uniform at higher reductions, at low 
reductions it is non-uniform and exhibits multiple peaks. It is 
interesting to note the peak point in the normal pressure curve 
at the entry which la the characteristic of all distributions. 

The normal pressure, as shown in figure 3.18, increases 
rapidly with decreasing values of R/h^^ and becomes less smooth. 
Again the interesting observation is at low values of R/h^^ there 
is a distinct peak at the entry which is not to be seen at high 
values of R/h^^ , for example, R/hj^ = 130. However, at the exit the 
variation is smooth and uniform for all values of R/h^^. 

Figure 3.19 shows the distribution of shear stress along the 

arc of contact from exit to the entry at different reductions. 

Since |t j=f|t I, at higher reductions t is more uniform similar 

to t^. From figure 3.20, where the distribution of shear stress 
n 

at different values of R/h^^ is shown, it is clear that the peak 
values of t^ also, similar to t^. Increase rapidly with 
decreasing R/h^. 

The variation of roll separating force with R/hj^ at 
different reductions is shown in figure 3.21. The figure suggests 
that while the roll separating force decreases with Increasing 
R/hj^, at all reductions in the range, this decrease becomes more 
rapid at increasing reductions. This is consistent with the trend 
of t^ with R/h^. This may be attributed to the indenting effect 
of roll into the work piece which increases with decreasing R/h^. 
At a given R/h^, roll separating force increases with reduction 
in spite of decrease in t^. So the increase in contact area seems 


50 



to be having more effect than decreaaing t^ on the roll 
separating force. Sioiilar observations can be made in the figure 
3.22 which shows the variation of roll torque with R/h^ at 
different reductions. However, in the case of roll torque the 
decrease with R/h. at higher reductions is more monotonous. The 
continuous increase in roll torque with reduction at any given 
R/hj^ can be justified by referring to the variation of shear 
stress in figure 3.19, where it is clear that the percentage 
length of the backward slip zone (distance from the entry to the 
neutral point ) at any given reduction increases with increasing 
reduction. Thus it follows that since the shear stress is 
positive in thebackward slip region roll torque increases with 
increasing reductions. 

3.2.2. Effect of the coefficient of friction 

The figure 3.23 shows the distribution of normal pressure 
along the arc of contact at different coefficients of friction. 
Uhile the shape of the curve and smoothness are restored for all 
values of coefficient of friction, the peak magnitude increases 
with increasing coefficient of friction and peak seems to shift 
towards the entry. Figure 3.24 which gives the variation of shear 
stress along the arc of contact at different coefficients of 
friction. It can be noted from this figure that the peak 
magnitude increases with increasing coefficient of friction and 
that the neutral point shifts towards entry. 

Figures 3.25 and 3.26 show the variation of roll separating 
force and roll torque with respect to the coefficient of 
friction, respectively. The modes of variation of roll force and 


51 


Ijrutv: 



roll torque In these graphs can be appreciated by referring to 
the figures 3.23 and 3.24 where the nomal pressure and shear 
stress distributions are depicted, respectively. It is 
interesting to see that since the contact area remains the same 
for all the three cases for which the curves are drawn, the 
variations, whatsoever, in roll force and roll torque are brought 
about by the change in shape and mode of variation of normal 
pressure and shear stress curves. From figure 3.25 it is evident 
that the roll force first increases with the coefficient of 
friction, however gradually, and then becomes uniform. This seems 
to be the consequence of the variation of normal pressure 
distributions which seem to become similar both^shape and 
magnitude with increasing coefficients of friction. Further the 
decreasing trend of roll torque curve with increasing coefficient 
of friction after reaching a peak can be appreciated by observing 
the figure 3.24. Here it is very clear to observe that with 
increasing coefficient of friction not only the neutral point 
shifts towards the middle of the arc of contact but the 
difference between negative and positive average values of shear 
stress decreases, both leading to a net decrease in the roll 
torque with Increasing coefficient of friction. However, 
meanwhile the fact that the computer program in this analysis 
ceases to converge for the values of coefficient of friction 
beyond 0.29 should also be pointed out. 

3.3. Effect of Strain Hardening 

In this section the effect of strain hardening is brought 
into light by selecting two sets of process variables. In the 
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firet aet the proceaa varlablea are auch that they are expected 
to give a minimum value of the deaign parameter to be obser ed 
and in the aecond aet they give an expected maximum of the aame. 
Theae expectationa are baaed on the foregoing parametric atudy. 
The deaign parametera are the aame which were conaidered for 
parametric atudy. Uhile the diacuaaion for the former two ia 
given by meana of grapha , the later two are preaented in the 
table 3.2. The percentage deviation for peak magnitudea of normal 
preaaure and ahear atreaa are alao included in the table 3.2. 

The figure 3.27 containa the diatribution of normal preaaure 
without and with atrain hardening for the two aeta of proceaa 
variablea. Uhile it la clear to aee that the normal preaaure la 
higher with atrain hardening than without atrain hardening all 
along the arc of contact in both the aeta, the peak deviation 
occura near exit in the firat aet whereas it occurs in the middle 
of the arc of contact in the aecond set. Further the percentage 
error in the estimation of t^ when the atrain hardening effect is 
neglected ia more for the aet 2 (see the table 3.2). It can be 
noted that the variation of t^ along the arc of contact ia more 
uniform for the aet 1 than for set 2. Figure 3.28 shows the 
effect of atrain hardening on the ahear atreaa diatribution along 
the arc of contact for the aame two seta of process variables. It 
shows that the percentage deviation ia more for aet 2. The 
underestimation of roll force and roll torque by not considering 
the effect of atrain hardening ia preaented in table 3.2. 
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3.4. The Study of Contours 

The figures 3.29 to 3.32 show the contours for the 
nor 0 ial~strain rate, equivalent strain~rate, deviatoric stress 
invariant and the equivalent strain in the domain for four sets 
of values of the three process variables. The four sets are so 
selected that the comparison with the contours of the first set 
of any of the successive three sets gives the effect of one 
process variable on the distribution of that parameter in the 
domain. The graphics routine for plotting these contours has been 
developed along with the computer program for the main finite 
element analysis. 

The figure 3.29 shows the contours for normal strain rate 
for steel for the four sets of process variables in Ca),(b),(c) 
and (d). The normal strain rates, at any given set of process 
variables first Increase from a small value as the metal enters 
the deformation zone and reach a maximum and thereafter decrease 
gradually as the metal exits from the deformation zone. However, 
a marked difference can be observed between these trends along 
the surface where the work piece is in contact with the roll and 
along the axis of symmetry. Uhile the normal strain rates along 
the axis of symmetry uniformly increase, reach a maximum and then 
decrease to a small value, just as explained above as a typical 
nature, along the work piece-roll interface they undergo a 
non-uniform variation with multiple maximums and minimums from 
inlet to the exit. This variation from uniform along the axis of 
symmetry to non-uniform along the interface la gradual. Because 
of this inhomogeinity of deformation pockets of normal strain 
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f fcB iTftdul't &loti^ And Jufict 1^6lov iTiiG in'tfii?£&cfi. It c&n 
observed from the figure 3.29 (a), (b), (c) and (d) that with 
increased reduction ratio ( a to b ) the homogeinity of 
deformation increases and with increased R/h^^ and coefficient of 
friction the inhomogeinity increases ( a to c and a to d, 
respectively ). The normal strain rates in 3.29Ca) first 
increase and reach a maximum value at the first contact point. 
Then it decreases towards the middle of the arc of contact 
reaching a minimum. It then again increases and attains another 
peak value and continuously decreases towards exit. They have two 
peaks on the surface and a minimum, Uith increased reduction, as 
shown in 3.29(b), the number of peaks reduce to one which occurs 
at a point away from the first contact point. Further the normal 
strain rates continuously increase till the peak is reached and 
then continuously decrease till they reach the cut off value at 
the exit. The increasing coefficient of friction and decreasing 
R/h^ seem to have little effect on the distribution of normal 
strain rates, as can be seen from 3.29(c) and 3.29(d). However, a 
minimum is introduced before a contour goes down to cut off value 
at the exit. Similar observations are possible from the 
equivalent strain rate contours shown in the figure 3.30. 

The contours of deviatoric stress invariant which is a 
measure of stress at any point required to cause a given amount 
of plastic strain at that point . In general, at any given set of 
process variables, the equivalent deviatoric stress first rapidly 
increases at the entry as the loading occurs, continues to 
increase into the deformation zone reaching a maximum somewhere 
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1ft the middle aftd theft gradually decreaaea towards exit. However, 
at the exit where unloading of the work piece occurs , the 
decrease In the devlatorlc stress invariant is again very rapid 
as the metal comes out of the deformation zone. Although the 
effective devlatorlc stress remains more or less uniform across 
the rolling direction in the work piece at the loading and 
unloading points the variation between and surface particles is 
distinct. This observation la common to all sets in the figure 
3.31 except that the variations at the inlet and exit are not 

shown in any of these figures. Comparison of figure 3.31Cb) with 
figure 3.31(a) suggests that the devlatorlc stresses become more 
or less uniform across the section with increased reductions. 
However, there is not much effect of other two process variables 
on the distribution of these contours. The figure 3.32 shows the 
contours of equivalent strain, also in four sets. The following 
observations can be made from the figures for the contours of 
equivalent strain. At the entry as the metal enters into to the 

deformation zone the strain increases rapidly from a very low 

value. In this zone clearly the surface particles undergo higher 
amounts of strain than those inside. As the metal progresses 
towards the exit, however, the inside particles catch up and 

trend becomes almost reverse. Uhile the trend at higher reduction 
ratios is more or less similar, because of longer path available 
for the particles to travel from the inlet to the exit this 
variation is gradual. Equivalent strains at higher coefficient of 
friction are shown in figure 3.32(a) keeping the other variables 
constant. This illustrates that at higher friction the surface 


S6 



&re too highly stf’fl.lned to be reached by the Ineide 
particles and therefore throughout the domain the equivalent 
strain across the sections increases from inside to the surface. 
This variation is gradual from a low value to a high value. 

3.5. Conclusions 

Based on the above discussions the following conclusions can 
be drawn. In the present analysis cold plane-strain flat rolling 
problem is solved by the finite element method at steady state by 
considering the material to be rigid-plastic strain hardening. 
The model can give good results over a range of the process 
variables viz. the reduction ratio, R/h^^ and the average 
coefficient of friction. The model has the capability of 
accommodating variable coefficient of friction which can be 
considered from the experimental results to obtain better 
results. In this work a detailed parametric study has been 
conducted to manifest the effect of the process variables on the 
important design parameters of rolling process. The error of the 
assumption of perfectly plastic material is highlighted. Contours 
have been shown to illustrate the distribution of stresses, 
strain rates and strains for an understanding of the physical 
process of rolling. 

3.7. Scope for Further Research 

As has been indicated in the above lines, the model is 
certainly amenable to the inclusion of the variable coefficient 
of friction which can further Improve the accuracy of results. It 
should also be pointed out that, a change in the mesh by bringing 
in smaller elements at appropriate places is likely to improve 
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results. Uiiile three dlotenslonal trestmejit is not sn 
improvement over every 2-D analysis( especially when the 
width/ thi ckness is very high the cost incurred in 3-D analysis 
can prove to be too high to appreciate the extra accuracy 
achieved), by carrying 3-D analysis some additional features like 
profile of the strip can be brought into light and the problems 
of shape rolling where the flow of metal in the third direction 
is significant to achieve the required shapes can be dealt with. 
The assumption that the metal is rigid at the entry and exit and 
that the roll is rigid may also be eliminated to predict better 
results . 
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Metal : Steel — Ana. 



Fig. 3-1 Comparison ot Analytical and Experimental 
Roll Force 


Metal'. Steel — Ana. 



Fig. 3-2 Comparison ot Analytical and Experimental 
Roll Force 
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Metal; Aluminium — Ana. 

t =0-10, R/h-j =12*5 ,h‘|=6-274 mm □ Exp. [^5] 


Reduction (7o) 

Fig. 3-3 Comparison ot Analytical and Experimental 
Roll Force 


Metal : Aluminium — Ana. 

t=0-10,R/h'i=39,h^=2-032mm □ Exp.[25] 
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Reduction (7e) 

Fig. 3-4 Comparison ot Analytical and Experimental 
Roll Force 






Fig. 3-6 Comparison of Analytical and Experimental 
Roll Force 
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Fig.3-8 Comparison ot Analytical and Experimental 
Roll Torque 
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Metal : Aluminium Ana. 

£ t = 0-10, R/h^ =12-5, h'|=6'276 mm ° Exp.[25] 



Fig.3-9 Comparison of Analytical and Experimental 
Roll Torque 


Metal Aluminium — Ana. 



Fig.3-10 Comparison of Analytical and Experimental 
Roll Torque 
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_ Metal '.Copper — Ana. 

E ^ "O'OG , R/h-i=12'5, hi=6-35mm □ Exp. [25] 



Fig.3'11 Comparison of Analytical and Experimental 
Roll Torque 



Fig.3'12 Comparison of Analytical and Experimental 
Roll Torque 
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^ Normal Pressure (GN/sq m) 

Normal Pressure (GN/sq.m) 


Metal : Aluminium __.Ana. 

t = 0-10,R/hi=12-5,h'i=6*274,7orr14*17— Exp.[25J 



•13 Comparison ot Analytical and Experimental 
Normal Pressure 


Metal- Aluminium — Ana. 

t = 0-10, R/h^ = 39,hf 2032, 7or=337 5 —Exp. [25] 



Fig. 3-14 Comparison ot Analytical and Experimental 
Normal Pressure 
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Arc oi conlQct (degrees) 


Fig.3-15 Comparison of Analytical and Experimental 
Normal Pressure 


Metal : Copper --Ana. 

1=0-06, R/hi=50,hiJ-562,%r=22-76 — Exp.[25] 



Fig.3-16 Comparison oi Analylicol and Experimental 
Normal Pressure 
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Metal : Steel 8 °/c 

f = 0.06, R/h|= 65, R = 65 mm 16°/ 
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Fig,3j7Distribution of normal pressure at different reductions 



R/hi =35 

Mela\ •• S^eel _____ =95 

f = 0.06, °/or = '6 R/hi= '3’® 
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(Exit) (Entry) 

Fig. 3.19 Distribution of shear stress at different reductions. 





Metal ; Steel 

f = 0.06, %r = 16, R= 65 mm 
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Metal ; Steel 
= 0.06 , R = 65 mm 
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R/h| 

Fig.5.21 Variation of roll separating force with R/hf at d 
reductions. 



Roll torque/Unit width{ kN-m/m) 


Metal : Steel f= 0.06, R= 65 mm. 



Fig. 3. 22 Comparision of variation of roll torque with R/h| at 
different reductions. 




Metal : Steel 

^hi = 65 , R = 65 mm 
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Fig.3.23 Distribution of normal pressure at different coefficients of friction 




Metal : Steel 
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Fig3.27 Distribution of normal pressure with and without strain hardening 
(a) minimum and (b) maximum in the range . 
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Fig. 3.28 Distribution of shear stress with and without strain hardening (a) minimum and 
(b) maximum in the range 





Tabl* 1. Material Propertiaa and Straaa-Strain Relationa 


Material 

Denaity 

In the Relation, 


■MM 


(kg/m^) 

C MPa ) 

n 

1/ K 

Reference 


Aluminium 

2800 

50.3 

0.260 

0.050 

Al-Salehi [ 25 

] 

St eel 

7870 

358.0 

0.300 

0.044 

Shida [ 26 ] 


Steel 

7870 

324.0 

0.295 

0.052 

Shida [ 26 ] 


Copper 

8600 

70.3 

0.490 

0.022 

Al-Salehi [ 25 

] 


Table 2. The Percentage undereatinat ion in the analyaia without 
atrain hardening for roll aeparatlng force, roll torque, peak 
magnitude of normal preaaure and abaolute peak magnitude of ahear 
atreaa 


Parameter 

Set 


Parameter Value 

% Error 

(w.r.t. the value 
with atrain hard.) 




1 

without 

with 

Roll Force 
CkN/m) 

Set 

1 

596 

707 

15.70 

Set 

2 

2075 

3111 

33.45 

Roll Torque 
CkN-m/m) 

Set 

1 

1.06 

1.15 

07.82 

Set 

2 

11.13 

16.35 

31.93 

(t ) 

^ n max 

( MPa ) 

Set 

1 

382 

485 

21.24 

Set 

2 

700- 

1160 

39.66 

I(| 

( MPa ) 

Set 

1 

24 

30 

20.00 

Set 

2 

135 

230 

41.30 
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^ no. of nodes 261 
no. of elaents 56 




(b) tr = 32 R/h^ = 130 R = <5 mm f = 0.06 



(c) tr = a R/hj =65 R = 65 mm f = 0.06 



Fig. 3.29 Normal Strain-Rate Contoure at Four Different Seta of 
Proceaa Variablea, (a), (b), (c) and (d). 
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Fie- 3-30 Equivalent Strain-Rate Contours at Four Different Sets 
of Process Variables, (a),(b),(c) and (d) 
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^ no. of nodes 261 
^ no. of elements 56 
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